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Abstract 

We compute the spectral function ^(k, uj) of a model two-dimensional high-temperature super- 
conductor, at both zero and finite temperatures T. The model consists of a two-dimensional BCS 
Hamiltonian with d-wave symmetry, which has a spatially varying, thermally fluctuating, complex 
gap A. Thermal fluctuations are governed by a Ginzburg-Landau free energy functional. We as- 
sume that an areal fraction cp of the superconductor has a large A (/3 regions), while the rest has 
a smaller A (a regions), both of which are randomly distributed in space. We find that A(k, cj) 
is most strongly affected by inhomogeneity near the point k = (vr, 0) (and the symmetry-related 
points). For c/j ~ 0.5, A(k, w) exhibits two double peaks (at positive and negative energy) near 
this k-point if the difference between Aq and A^ is sufficiently large in comparison to the hopping 
integral; otherwise, it has only two broadened single peaks. The strength of the inhomogeneity 
required to produce a split spectral function peak suggests that inhomogeneity is unlikely to be 
the cause of a second branch in the dispersion relation, such as has been reported in underdoped 
LSCO. Thermal fluctuations also affect j4(k, w) most strongly near k = (vr, 0). Typically, peaks 
that are sharp at T = become reduced in height, broadened, and shifted toward lower energies 
with increasing T; the spectral weight near k = (vr, 0) becomes substantial at zero energy for T 
greater than the phase-ordering temperature. 
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I. INTRODUCTION 



In the last several years, the measured electronic properties of cuprate superconductors 
have shown considerable evidence of inhomogeneities. For example, spatial variations of the 
superconducting energy gap and of the local density of states spectrum have been observed 
in scanning tunneling microscopy (STM) experiments l|, [3, I^, 0, [s, 6, 7, ^. There have also 
been a number of reports of magnetic and charge ordering in these materials, which also 



indicate inhomogeneities [9|,ll0|, 
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la, llTJ • Other studies of cuprates have 



shown that electronic states within certain energy ranges show checkerboard-like spatial 
modulations [isl - 



191, 



A number of theoretical 

d, y, t 
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approaches have been developed to model these inhomogeneities 



extended in a recent article 
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33( 1 ■ These works are reviewed and 



34 to 



3J]. In the present article, we use the approach of Ref. 
explore how the spectral function of a d-wave superconductor is affected by gap inhomo- 
geneities and thermal fluctuations. 

The spectral function of cuprate superconductors has been studied theoretically by a 
number of groups, though most have omitted the effects of quenched inhomogeneities. For 
example, Wakabayashi et al. 35|] used a weak-coupling BCS theory combined with a Green's 
function approach to explain the narrow quasiparticle peak at the gap edge, which has been 
observed by ARPES experiments in overdoped cuprates along the antinodal direction. Fieri 
et al. Issj] , using a Nambu formalism, have studied a model which includes pairing fluctuation 
effects, and which accounts for some features of the single-particle spectral function as 

n 

observed in certain cuprates. Zacher et a/.[37[ have used a cluster perturbation technique 
to compute the single-particle spectral function of the t-J and Hubbard models, and to 
study stripe phases in the cuprates. Faramekanti et al. 38|] have studied a Hubbard model 
for (uniform) projected ci-wave states. They used a variational Monte Carlo technique in 
which one of the variational parameters is the magnitude A of the pairing field, and find 
that A, as a function of doping, scales with the (vr, 0) hump and T* as observed in ARFES 
experiments. 

In another recent work, using a generalization of BCS theory, Chen et al. 39|] found a 
sharpening of the peaks in the spectral function as T is reduced below T^, similar to what 
we find as discussed below. (Their model involves a homogeneous superconductor.) Hotta 
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et al. 



40[ | have used a self-consistent t-matrix approximation to study a model for s and d 



wave superconductivity at finite T. They found a gap in both the single-particle density of 
states and the spectral function even above the superconducting transition temperature T^, 
the energy scale for the pseudogap is found to be the Cooper-pair binding energy. 

More recently, Mayr et al. 4l|] have introduced an extended Hubbard model, which in- 
cludes both superconductivity and antiferromagnetism; they found that quenched disorder 
is a necessary ingredient for that model to reproduce the double branch, or split band, ob- 
served in angle-resolved photoemission (ARPES) experiments on La2_2:Srj.Cu04. Finally, a 
model to study how A(k, u) is affected by thermal fluctuations of the phase, but not the 
amplitude, of the superconducting order parameter has been treated by Eckl et al. |42|] for 
homogenous systems. 

In the present work, we propose a simple model for A{k, uj). This model can exhibit a split 
peak near k = (vr, 0), but only for certain parameter choices which are unlikely to be realized 
experimentally. Our model consists of a BCS superconductor with (i-wave symmetry, where 
the pairing field (given by the superconducting order parameter) is inhomogeneous, and 
is also subject to thermal amplitude and phase fluctuations at finite T. We assume that 
those thermal fluctuations are governed by a discretized Ginzburg-Landau (GL) free energy 
functional. To compute the spectral function, we use exact numerical diagonalization of the 
BCS Hamiltonian on a finite lattice and average over many different configurations of the 
thermally fluctuating superconducting order parameter as obtained using the Monte Carlo 
technique. Thermal averages are obtained by averaging over these configurations. 

Inhomogeneities are introduced in our model phenomenologically. The atomic lattice is 
subdivided into cells, which we call XY cells, of size 2x2 atomic sites; within each such cell, 
we assume A to be constant. Then we choose the coefficients of the GL free energy functional 
so as to give, at T = 0, a binary distribution of the superconducting order parameter |A| at 
each atomic lattice site. XY cells with small and large |A| values are called a and (3 cells, 
respectively. We take the distribution of a and (3 cells on the atomic lattice to be random, 
as suggested by STM experiments. The two parameters which we vary in our calculations 
are (i) the area fraction C/3 of (3 cells, and (ii) the magnitude |A^| of the gap in those cells. 
The value of |Aa| is kept the same through our calculation, and is inferred from the STM 
experiments. 

At T = 0, we find that the main consequence of this binary, random distribution of A 
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is to broaden the peaks in the spectral function A(k, u) near the points k = (vr, 0) [and the 
symmetry-related points at k = (— vr, 0) and (0, ±7r)]. This broadening is most pronounced 
at Cj3 = 0.5. For a sufficiently large ratio A/j/A^ of the large to the small gap, we find that 
y4(k, cu) near k = (vr, 0) shows two peaks rather than one ("split band regime"). Otherwise, 
we find a single peak which is broadened by disorder. Although the experimental ARPES 
results of Yoshida et al. 43| also show a double peak, there are several reasons to believe 



that this split peak is not caused by the kind of inhomogeneities we consider here. This 
point is discussed further below. 

At finite T, we find that, near k = (vr, 0), the originally sharp coherence peaks of A(k, ti;) 
as a function of uj broaden and shift to lower energies with increasing T. This broadening is 



similar to that found in the calculations of Eckl et al. [42j, which omits quenched disorder 
and also include thermal fiuctuations only in the phase but not the amplitude of the gap. 
These calculations focused on T near that of the phase ordering transition. By contrast, we 
present calculations showing how y4(k, cu) evolves near k = (tt, 0) as a function of uj over a 
broad range of temperature, including both amplitude fiuctuations and quenched disorder. 
The remainder of this paper is organized as follows. In Section II, we briefiy describe our 



model, which is already presented in Ref. 



34 . In Section III, we give our numerical results. 



followed by a discussion and conclusions in Section IV. 



II. MODEL 



A. Microscopic Hamiltonian 

We consider the following Hamiltonian: 

H = 2^^ Ujclcj^ + 2 ^{AijdiCj^ + c.c.) - f^Yl ^^^'^^^'^ 

Here, denotes a sum over distinct pairs of nearest neighbors on a square lattice with 

sites, c]^ creates an electron with spin a (| or I) at site j, /i is the chemical potential, Ajj 
denotes the strength of the pairing interaction between sites i and j, and tij is the hopping 
energy, which we write as 

tij thop- (2) 

where thop > 0. 

4 



Following a similar approach to that of Ref. 3J] and Ref. 4J] we take to be given by 

I A, I 



A,; 



1 |A,| 
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where 



+ 6j)/2, if bond (z, j) is in x-direction, 
+ 9j)/2 + vr, if bond (i, j) is in y-direction, 



and 



A. 



|A,|e 



(3) 



(4) 



(5) 



is the value of the complex superconducting order parameter at site j. The sums in ([T]) are 
carried out over a lattice we will refer to as the atomic lattice (as distinguished from the 
lattice, described below). The first term in eq. ([1]) corresponds to the kinetic energy, 
the second term is a BCS type of pairing interaction with (i-wave symmetry, and the third 
term is the energy associated with the chemical potential. 



B. Numerical Calculation of Spectral Function 



We wish to compute the spectral function A(k, uj) for the system described by the Hamil- 
tonian ([T]). Given the Aj's, tjj, and /x, y4(k, u;) is computed through 



A(^,k,{A,})= [\Unm^5{uJ~En) + \Vn{\^)\''5{uJ + En)i 



n,E„>0 



where 



i=l 
1 ^ 



En is the nth eigenenergy of Hamiltonian ([T]), and 



Un[ri) 



is its nth eigenvector, as described in detail in 



3^. Here, 



(6) 

(7) 
(8) 



(9) 



(10) 
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and Hi and are integers in the range [0, A^^. — 1] and [0, A^y ~ !]• our numerical 
calculations, we take the size of the atomic lattice to be = N^Ny, where x and y are unit 
vectors in the x and y directions, and is the lattice constant. We use periodic boundary 
conditions, ^/'„(r) = ipn{r + N^aox) and V'n(r) = ^/'„(r + Nyaoy), which leads to k- vectors of 
the form 

9-n" / m m \ 

(11) 



27r / rrix rriy 
k = — —X + -j-y 



34 Basically, we subdivide 



The detailed procedure to obtain Aj is described in Ref. 
the atomic lattice into cells, which we call XY cells, of size x 'Co- Here .^o is the T = 
Ginzburg-Landau (GL) coherence length, which we take to be an integer multiple of a^. The 
value of Aj is assumed to be the same for each atomic site within a given XY cell, and is 
governed by the following discretized GL free energy functional: 



- = E( — -1 



M 



i=l 



T 



A, 



M 

- 

18.76 A2( 



A, 



A, 



A, 



cOi 



Xj{0)kBTcOj 
(12) 

Here Ki = ^'^(i/[32(9.38)7rm*'^/i^], where m* = 2me is twice the electron mass, /is is the 
Bohr magneton, and d is the thickness of the superconducting layer. If d = lOA, Ki = 2866 
eV A^. Aj is the complex gap parameter in the i*'* XY cell. In eq. f ll2l) . the sums run over 
the lattice of XY cells, each of which contains {^/clqY atomic sites. 

We choose the coefficients of this GL free energy functional Tcoi and Ai(0) to have binary 
distribution on the XY lattice, corresponding to either a small or a large value of |Aj|. We 
call an XY cell with a small (large) value of |Aj| an a {(3) cell, while the area fraction of j3 
cells is called c^j. The corresponding values of Tcoi and Aj(0) are denoted TcOa and Tc^p. At 
T = 0, in a homogeneous system made up entirely of a {(3) cells, the magnitude of |Aj| will 
be the same in each XY cell and given by the minimum of the corresponding free energy 
functional F, i. e. |Aj| = v^Q^fceT'cOa ("s/O^fc^Tco/j). In the binary case (0 < c^j < 1), at 



T = 0, we will still generally have |Aj| = vQ^/^^TcOi, although this value may be modified 
slightly by the proximity effect term in F / Ki [the last term in eq. f|T2l) ]. 

We compute y4(k, cj) at T = by diagonalizing the model hamiltonian (1) using Aj 
determined by minimizing the Ginzburg-Landau free energy F. This minimum value will 
always correspond to gaps A,, = |Aj|e*^' such that all the phases 6i are equal. At finite T, we 
compute A(k, cu) as an average over different configurations {A^}. These are obtained, as 



in Ref. 



34l |. by assuming that the thermal fiuctations of the Aj are governed by the GL free 



energy functional F described above. Thus, F is treated as an effective classical Hamiltonian 
and thermal averages such as (y4(k, tu)) are computed as 

(^(k, -)) = Jnl.d^^^e-^I^.T • (13) 

We will be using the GL free energy functional at both T = and finite T in spite of the 
fact that it was originally intended for T near the mean-field transition temperature. Strictly 
speaking, the correct free energy functional near T = should not have the GL form but 
would be expected to contain additional terms, such as higher powers of We use the 
GL form for convenience, and because we expect that it will exhibit the qualitative behavior 
that would be seen in a more accurate functional - that is, the effects of inhomogeneities 
would be quantitatively the same in the GL model as in a more accurate model containing 
additional powers of IV'P- 

To obtain A(k, uj) for a given distribution of the Aj's, we diagonalize the Hamiltonian (1) 
for that configuration, then obtain A(k, cu) using eq. (6). The canonical averages are then 
evaluated using a Metropolis Monte Carlo technique to determine the canonical distribution 
of the Aj's at the temperature of interest. The detailed description of this Monte Carlo 



approach are given in Section IV. A of Ref. 3J]. As noted there, we first choose the values 



of TcOi and Aj(0) in each XY cell, taking these to be quenched variables. In contrast to our 
calculations of Ref. sj] , we do not include a smoothing magnetic field to reduce finite-size 
effects; as a result, our results have more numerical noise than do our earlier results for the 
density of states. 

For the present model calculation, we arbitrarily set the chemical potential /i = 0, for sim- 
plicity, and use thop as the unit of energy. This corresponds to half filling in the band model. 



Exactly half filling would correspond to x = in Lai.Sri_^Cu04 (LSCO), for example |45[|. 
It should be noted that some of the most interesting experimental results for the spectral 
function I43I are carried out in the underdoped superconducting regime of the phase diagram, 
where /i is slightly negative. If we set yU 7^ in our model, this leads to unequal integrated 
weights of the spectral function peaks at positive and negative energy, but we have found 
that otherwise our numerical results are not very different from those at /i = 0, for our 
model Hamiltonian. However, the present results and model, for reasons which we discuss 
below, are probably not directly relevant to those experiments. 

In order to show that our results are not strongly affected by setting = 0, we have also 



done simulations using /i 7^ 0. For example, in Figure [8] we present results using fi = —0.05. 
For this value of /i, the average number of electrons per site, defined as, 

1 ^ 

i^) = NY.(^^) (14) 

i=0 

with 

ni = ^4^Q,<, (15) 
is found to be (n) ~ 0.94. This corresponds to a strongly underdoped cuprate x ~ 0.06. 

C. Homogeneous systems 

For a homogeneous system at T = 0, Aj = A and we can rewrite Hamiltonian ([1]) as 

H = Y1 ^kicka + ^(AkCk|C_kT + c.c.) - ^ ^ icka, (16) 

k,(7 k k 

where ek = —2t[cos{kxao) + cos(fcyao)] and Ak = ^A[cos{kxao) — cos(A;yao)]. In obtaining 
(fT6|) we have used = ^,^exp(— ik ■ r^) and its hermitian conjugate. In this case, 
the excitation energies of the system are given by [47] 

= 7(6k - ^)2 + A2 (17) 

The corresponding spectral function will be a sum of two delta functions, as indicated by 
eq. 

III. NUMERICAL RESULTS: INHOMOGENEITIES AND THERMAL FLUCTU- 
ATIONS 

In this section we present our numerical results for A(k, u) for inhomogeneous systems 
both at zero and finite temperatures; for reference, we also show the corresponding results 
for homogeneous systems in some cases. For T = 0, we use 48 x 48 atomic lattices used, 
while at finite T we used lattices of 32 x 32. In al cases, we use 2x2 XY cells. Through the 
rest of this article, we show energy measured in units of thop, distance in units of ao, and k 
in units of l/flo- 
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A. Zero temperature 

Before describing our results at zero temperature, we first comment on our choice of 
gap parameters used in the calculations. Our primary goal is to ascertain what kinds of 
qualitative spectral functions could result from the type of inhomogeneity described by our 
models, not to compare directly to experiment. For this reason, we will examine gaps which 
are, in general, substantially larger (in units of t/iop) than those which would describe realistic 
cuprate superconductors. This point is examined further in the discussion section. 

With this preamble, we now present our results at T = 0. Fig. [1] shows the spectral 
function A(k, u) (represented as a contour plot) as well as plots of the dispersion relation 
Eif_ as a function of k, for two homogeneous systems: one with A = 0, and another with 
A = 0.42. For such homogeneous systems, A(k, uj) is simply proportional to the sum of two 
delta functions: A(k, u) oc 6{uj—Ei^)+6{u!+Ei^). In parts (a) and (b), the dark (light) regions 
correspond to regions where E]^, as calculated from Eq. f|T7|) . is large (small); these are shown 
for all k vectors in the first Brillouin zone (BZ). For a system with A = [Fig. W^a)] there 
are four lines (white) in k-space for which E]f_ = 0: ky = ±/c^ ± vr. When A > [Fig. [Il^b)], 
the lines are reduced to four points: {k^ = Tr/2, ky = ±7r/2) and {k^ = — vr/2, ky = ±7r/2), 
located at the center of the white blobs in Fig. [TJ^b). In Fig. Wic) and (d), density plots of 
A(k, uj) as a function of u are presented for those homogeneous systems at selected k values. 
These values lie along three standard lines in the first BZ: from k = (0,0) to k = (vr,0), 
from k = (vr,0) to k = {tt,tt), and from k = (vr, vr) to k = (0,0). Dark (light) regions 
correspond to large (small) values of the spectral function. For each k in these homogeneous 
systems, there is a sharp peak in A(k, uj), whose energy and width are indicated as the very 
short dashed lines in the plot. Also, the spectral function is clearly most strongly affected 
by a finite value of A near k = (tt, 0), where an energy gap of magnitude A opens around 
u — fi = 0. 

Figures [2] and 3 show the spectral function of several inhomogeneous systems with differ- 
ent concentrations C/j of f3 cells, at T = 0. In these systems, the atomic cells within the f3 
cells have A = 1.26, and are randomly distributed in the atomic lattice, while a cells, which 
occupy the rest of the lattice, have AO. 42. Fig. H] shows a representative arrangements of a 
and (3 cells for an 16 x 16 XY lattice with = 0.1. We can observe in Fig. [2] that the dis- 
order introduced by this binary distribution of the superconducting order parameter affects 
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mostly the region k = (vr, 0). This disorder effect is almost unobservable for cp = 0.9: the 
results are almost the same as those for = 1.0, a homogeneous system with only j3 cells. 
On the other hand, a small but noticeable disorder effect is observed C/3 = 0.1, in the form 
of a slight broadening of the spectral function at k ~ (0.87r,0). However, it is the system 
with = 0.5 the one that shows a most dramatic blurring of the energy in the region of 
k = (vr, 0), as we now discuss. 

Since the effects of the binary distribution of A are more pronounced near k = (vr, 0), we 
have also plotted A(k, u) versus u for fixed k = (vr, 0), at different values of C/3 in Fig. [31 For 
the pure a system, C/? = 0.0, two sharp peaks appear at \uj\ = = 0.42. When a fraction 
0.1 of the a XY cells are replaced by (3 cells, C/3 = 0.1, the height of the peaks decreases from 
about 45 (arbitrary units) to about 15, with a corresponding broadening of the peak and a 
shifting of the weight toward a higher energy. At C/j = 0.5, the peak height is only about 1.5, 
and the width is very large; the peak fills the entire frequency range from u = Aq, = 0.42 
to A^ = 1.26. At Cj3 = 0.9, most of the weight of A(k, cj) shifts to u = A/3 = 1.26, with 
a slight broadening near the bottom of the peaks, which is, however, less pronounced than 
the corresponding broadening of the = 0.1 peaks. At = 1.0 the peaks become sharp 
at 00 = ±Ap = ±1.26. 

We now discuss the calculated effects of a binary gap distribution on systems similar 
to those of Fig. [2], but with A^ = 2.52 instead of A^ = 1.26. These systems, like the one 
previously discussed, have A^ = 0.42. For concentrations = 0.1 and cp = 0.9, the effect of 
inhomogeneities qualitatively resembles that seen for A/? = 1.26: they produce broadening 
of the spectral function near k = (7r,0). The main difference is that the broadening is 
slightly greater for A^ = 2.52. However, the case cp = 0.5 shows a real qualitative change: 
the spectral function now splits into two well-defined peaks for k vectors near (vr, 0). We 
can better visualize this effect by looking at Fig. [6], where we plot A(k, u) versus u for fixed 
k = (vr, 0) and several values of C/3. Clearly, A(k, cj) for C/3 = 0.1 and C/j = 0.9 behaves 
similarly to the case A^ = 1.26: slightly broadened peaks at an energy near the A of the 
majority of the XY cells, i. e., at u = 0.42 for c/3 = 0.1 and aX u = 2.52 for C/3 = 0.9. But for 
C/3 = 0.5, A(k, cij) shows several peaks, two of which are particularly clear: one at ~ 0.42 
and the other at u ~ 2.52. This is the "split band" regime one expects for large contrast 
between A^^ and A/3. 

In order to better visualize how the spectral function depends on disorder, we have 
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calculated A(k,u!) as function of \Aa\ for a fixed ratio \A[^/Aa\ = 6, at C/3 = 0.5. The 
results are shown in Figs. [7] and 8. This series of plots clearly shows the evolution of 74(k, cu) 
from a split-band regime at |A^| = 2.52 (in units of thop) to a broadened single band for 
|A^| = 0.22\thop) or smaller. In general, we find that the split band regime occurs only if 
the difference | — |Aq| is of order thop or larger; otherwise, A(k, u;) at k = (vr, 0) is simply 
the sum of two broadened peaks at positive and negative energies. 



B. Finite temperatures 

Fig. [9] shows the T-dependence of A(k,uj) for a system with c/3 = 0.1. The a and /5 cells 
are now characterized by values of tco such that at low T, A = 0.42 in the a cells, and 
1.26 in the jS cells. The value of A itself at finite T will, of course, thermally fluctuate, as 
governed by the GL free energy functional discussed at the end of Section II. The spectral 
function presented here is therefore an average of A{uj, k, {Aj}) over different configurations 



{Aj} obtained by a Monte Carlo sampling procedure, as described above and in Ref. 34 1. 
Hereafter, we denote this ensemble average simply as A(k, cu). 

Fig. [9] shows that, as in the case of quenched disorder, A{k,uj) is most strongly affected 
by thermal fluctuations near k = (vr, 0), where it broadens more and more with increasing T. 
In addition to this broadening, the peaks can be seen to shift towards smaller energies. This 
behavior can be seen more clearly in Fig. [101 which shows A(k, tu) at k = (vr, 0) as a function 
of u — fi. We observe that at T = 0, A(k, u) shows relatively sharp peaks at u — fi = ±0.42, 
with some disorder-induced broadening only in the wings of the peak. At t = 0.01, the peak 
height of A(k, u) decreases from ~ 17 to about ~ 7, with a correspondingly increased width. 

As the temperature is increased, the system eventually undergoes a phase-disordering 
transition, above which the superconductor loses phase coherence. For the parameters used 
in Fig.Hni this transition occurs at tc — 0.035. tc is the phase ordering transition temperature 
in units of thop/ks- We use a dimensionless temperature t = kBT/thop in these plots. At 
t = 0.03, near but slightly below the phase ordering temperature — 0.035, the height of 
the peak is further decreased, its width further increased, and its energy shifted to a still 
lower energy. At t = 0.05 > tc, the peak shifts still further toward lower energy, but the 
maximum remains at finite energy. 
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IV. DISCUSSION 



We have presented a simple model to study how the spectral function of a model d-wave 
superconductor is affected by quenched inhomogeneities and by thermal fluctuations of the 
superconducting order parameter. The model consists of a BCS Hamiltonian for an order 
parameter with (ia;2_j^2-wave symmetry, which has a position-dependent pairing field, and 
which also undergoes finite-temperature thermal fluctuations. The spatial dependence we 
assume for the pairing fleld is motivated by recent STM experiments on Bi2212: we assume 
two types of regions: and a region with a small gap, and a /? region with a large gap. 
To treat thermal fluctuations (of both amplitude and phase of the superconducting order 
parameter), we assume that they are governed by a suitable Ginzburg-Landau free energy 
functions, which we treat by classical Monte Carlo simulations. 

At T = 0, we flnd that A(k, cu) is most strongly affected by disorder near k = (vr,0). In 
general, this effect consists of a broadening of the peaks of A(k, lS) (plotted as a function 
of io for flxed k). However, at area fraction = 0.5, we flnd that quenched disorder can 
have two qualitatively different effects, depending on the relative magnitudes of Aq, and Ap. 
If the difference between Aq, and A^ is small, A(k, cu) has a single, broad peak for k near 
(7r,0), extending from ~ Aq to ~ A^. But for a large enough difference between Aq and 
A^, the A{k,uj) show a characteristic "split-band" behavior: instead of a wide, single peak, 
there are two prominent peaks, at u = A^ and u = A/j. 

Thermal fluctuations of the pairing fleld also have their strongest effect on A{k,uj) near 
k = (vr,0). The effect consists of a gradual broadening of the T = peaks with increasing 
temperature, and also a shifting of those peaks towards lower energies. However, no dramatic 
change is noticeable near the phase-ordering transition. 

Finally, we comment on the possible connection, if any, between our results and exper- 
iment. In recent angle-resolved photoemission studies by Yoshida et a/|43|, for LSCO, it 
was observed that for doping level x = 0.03, a second branch developed in the dispersion 
relation near k = (vr, 0). An explanation for the presence of these two branches has recently 
been suggested by Mayr et al. [41,]. These authors showed that the extra branch could be 
explained by a model with quenched disorder, in which the material breaks up into spatially 
separated superconducting and antiferromagnetic patches. 

In the present work, we flnd that a similar effect, with two spectral peaks, can be produced 
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if there are spatially distinct superconducting regions with sufficiently different supercon- 
ducting gaps. However, we also find that a split spectral peak can be produced only if the 
magnitudes of the gaps |Aq| and |A^|, and of their difference, is much larger than seems 
physically reasonable. Specifically, unless one of the gaps in the bimodal distribution is 
around 2.5thop, we do not obtain a split peak in A(k, u) at the point k = (tt, 0). For typical 
values (thop ~ 200 meV), this would represent a |A| of around 0.5eV. Since the average 
value of |A| in most of the cuprate superconductors is ~ O.OSeV, it seems most unlikely that 
random spatial fluctuations in |A|, due to quenched disorder, could produce such a large 
gap locally. Furthermore, even with such large quenched fluctuations in the gap, we need a 
bimodal gap distribution to obtain a split spectral function - equally large quenched fluctua- 
tions, but with a continuous distribution due to quenched disorder, would probably not give 
rise to a split spectral function. Therefore, it seems very improbable that our model could 
account for the second branch in the dispersion relation reported in Ref . [4^ . However, our 
results should give a reasonable picture of how quenched gap inhomogeneities affect A(k, u) 
in a d-wave superconductor over a range of parameters. 
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FIG. 1: Contour plots of the energy E'k of single-particle excitations, and the corresponding spectral 
function spectral function A(k, a;), for two homogeneous systems at zero temperature: one with 
A = 0, and another with A = 0.42. In parts (a) and (b), the dark (light) regions correspond to 
large (small) values of the excitation energies, as calculated from Eq. (jl7p . In parts (c) and (d), 
we plot the positions of the peaks of A(k, w) for these two homogeneous systems. 
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FIG. 2: Same as Fig. 1(c) and 1(d), but for several inhomogeneous systems, with different concen- 
trations of (3 (large-gap) cells, at T = 0. In all four plots, atomic sites within the a and /? cells have 
A = 0.42 and A = 1.26, respectively; the cells are randomly distributed over the atomic lattice, 
as illustrated in Fig. [J] for the case cp = 0.1. The dark (light) regions correspond to large (small) 
values of A(k, w). A more quantitative view of A(k, w) is shown in Fig. [3] for k = (7r,0). 
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FIG. 3: Spectral function A[k = (tt, 0), u] as a function of oo of systems with different concentrations 
cp of P cells at zero temperature. We have taken Aq, = 0A2thop, = 1.26thop 
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FIG. 4: A typical realization of disorder in a system with a concentration = 0.1 of /3 cells 
(white) immersed in a background of a cells (grey). Each cell contains four atomic sites. 
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FIG. 5: Same as Fig. O except that /3 cells have a A/j = 2.52 instead of A/3 = 1.26. 
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FIG. 6: Same as Fig. [3] but for /3 cells which have A = 2.52 instead of A = 1.26. 
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k=(7i,0), Ae=6A„, t=0.0, Cg=0.5 




FIG. 7: Spectral function A[k = (tt, 0),a;] as a function of ui for a system with concentration 
cp = 0.5 of /? cells at r = 0, plotted as a function of the magnitude |Aq,| of the component with 
the smaller gap. In all cases, |A^/Aq,| = 6. 
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k=(7i,0), Ap=6A„, t=0.0, Cp=0.5 




FIG. 8: Same as Fig. [71 but with chemical potential = —0.05 instead of /i = 0. With our model, 
setting fj, = —0.05 results in having an average occupation number (n) ~ 0.94, which corresponds 
to an strongly underdoped cuprate x ^ 0.06. 
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FIG. 9: Plots of the spectral function A ( k, u) for a system at zero and finite temperatures. The 
system has a concentration C/3 = 0.1 of /? cells with = 1.26 randomly distributed in a background 
of a cells having = 0.42. The phase-ordering temperature is tc ~ 0.035. The detailed evolution 
of the curve A[k = (vr, 0), w] versus lj can be seen in Fig. 
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FIG. 10: Spectral function A[k = (tt, 0),a;] as a function of oo for a system with = 0.1 at 
different temperatures. The phase-ordering temperature of the system is tc — 0.035. 
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